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Abstract 

For studying the group theoretical classification of the solutions of the 
density functional theory in relativistic framework, we propose quantum 
electrodynamical density- matrix functional theory ( QED-DMFT ). QED- 
DMFT gives the energy as a functional of a local one-body 4x4 matrix 
Q(x) = —(ip{x)ij){x)), where ip and ij) are 4-component Dirac field and its 
Dirac conjugate, respectively. We examine some characters of QED-DMFT. 
After these preparations, by using Q(x), we classify the solutions of QED- 
DMFT under 0(3) rotation, time reversal and spatial inversion. The behavior 
of Q{x) under nonrelativistic and ultrarelativistic limits are also presented. 
Finally, we give plans for several extensions and applications of QED-DMFT. 



* Corresponding author; Tadafumi Ohsaku, tadafumi@fuji.phys.wani.osaka-u.ac.jp 



1 



keywords; density functional theory QED, QED-DMFT, QED-TDDMFT, QED-CAS- 
DMFT, group theory symmetry and broken symmetry heavy elements. 



I. INTRODUCTION 



The density matrix was first introduced by Dirac [1], and after it was generalized by 
Husimi [2] and Lowdin [3] ( the reduced density matrix ) . This was recognized as a very useful 
tool to investigate interacting many-body theory and it gave various important concepts 
until now. The definition of the one-particle density matrix is given as 



Here -0 and i\y are 2-component electron field and its hermitian conjugate, respectively. 
On the other hand, the density functional theory ( DFT ) was first given by Hohenberg- 
Kohn-Sham [4], and it was developed by many researchers. People applied DFT to various 
many-fermion systems, not only many-electron systems like atoms, molecules and solids, but 
also nuclear systems. Today, DFT is well accepted, and used to investigate various physical 
and chemical characters of systems [5~7]. 

The symmetry aspect of DFT was also discussed by using the one-particle local 2x2 
matrix density: 



By using (2), Gorling examined the symmetry of the Hohenberg-Kohn theorem in de- 
tail [8]. Weiner and Trickey investigated the group theoretical classification of nonrela- 
tivistic DFT solutions [9], based on the work of Fukutome in which the classification for 
generalized Hartree-Fock ( GHF ) theory was presented [10]. Yamanaka et al. [11] ap- 
plied the generalized spin orbital ( GSO ) DFT to various molecular systems and studied 
the symmetry breaking phenomena in DFT solutions. We should mention also, it was 
demonstrated that the Hohenberg-Kohn theorem varid for one-particle density matrix func- 
tional E[y(x,x')] [12~15] . DFT for time- dependent phenomena was also investigated ( the 



~t{x,x') af} = (ipl(x')ip a (x)). 



(1) 



p(x) a/3 = (lpl(x)^ a (x)) . 



(2) 
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time-dependent density functional theory ( TDDFT ) ). Runge and Gross examined the 
Hohenberg-Kohn theorem for time- dependent external potential ( the Runge-Gross theorem 
) [16]. TDDFT was applied to estimation of electronic excitation or susceptibilities [17]. 

By the way, the first generalization of relativistic Hohenberg-Kohn theorem was intro- 
duced by Raj agopal- Callaway [18], and several attempts were made to investigate relativistic 
density functional theory ( RDFT ) [19~24]. Based on these studies, an important review 
of RDFT was presented by Eschrig et al. [25]. Today, the Hohenberg-Kohn-Sham scheme is 
well established for RDFT [26~28]. Relativistic time-dependent density functional theory ( 
RTDDFT ) was also investigated [29,30]. 

In this note we investigate quantum electrodynamical density-matrix functional theory 
( QED-DMFT ), and its group theoretical aspects. Recently, we constructed the theory of 
quantum electrodynamical self-consistent fields ( QED-SCF ) [31]. We derived the time- 
dependent Hartree-Fock ( QED-TDHF ) theory, Hartree-Fock condition and random phase 
approximation ( QED-RPA ). We also presented the group theoretical classification of quan- 
tum electrodynamical Hartree-Fock ( QED-HF ) solutions [32]. Based on our previous re- 
sults, in this work we apply the group-theoretical method to QED-DMFT. Our theory can 
become a basis or a starting point to investigate the phenomena of symmetry breaking in 
many-body system [9~11,33] under strong relativistic effects, by using DFT framework. In 
the work of Yamanaka et al. [11], the results indicate that the GSO treatment of DFT is very 
important and usuful for description of some polyradical clusters. Therefore, it is interesting 
for us to extend their treatment to the relativistic framework, and study the phenomena of 
symmetry breaking in heavy elements, analogous with the nonrelativistic GHF theory [10]. 

This paper is organized as follows. In Sec. II, we introduce QED-DMFT. Here, "matrix" 
means a matrix in spinor space. In fact, we do not use a non-local one-body matrix like 
7(2;, x') a p in nonrelativistic case, rather we intend to discuss a local one-body matrix like 
p{x)a/3 in nonrelativistic case for obtaining simple classification under group operations. 
The relation between various available DFT and QED-DMFT is also discussed. Quantum 
electrodynamical time-dependent density-matrix functional theory ( QED-TDDMFT ) is 
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also investigated. In Sec. Ill, the group theoretical classification of QED-DMFT solutions 
is discussed. The results can be applied to elucidate the phenomena of broken symmetry, 
especially in atomic systems under relativistic effects. Finally in Sec. IV, possible extensions 
and applications of QED-DMFT are discussed. Conclusion with summary of this work is 
also presented. 



II. QED-DMFT 



A. Introduction for QED-DMFT 



In this subsection, we introduce QED-DMFT. First, we give a brief summary of QED. 
We start from the following Lagrangian: 



C(x) = --F^(x)F^(x) + iP(x)(i^ - mo)ip(x) + e^(x)^(x)(A^x) + A e *\x)). (3) 

Here, is a ( virtual ) photon field which describes the interaction between fermions, and 
A e * 1 is a classical external field. F fll/ (x) = d fl A u (x) — dyA^x). m is the bare mass of 
fermion, ip(x) is the Dirac field, ifj(x) is its Dirac conjugate, x = (xo,x) indicates time-space 
coordinate. The 7-matrices are defined by ( under the standard representation ) 



<i 0^ 



7 



-/ 



/ ^ 



-o l 



(4) 



where I is the unit matrix, and Oi ( i=l,2,3 ) is the Pauli matrix: 



v i 0, 



—1 



i 



/ i ^ 



-1 



(5) 



The Hamiltonian for a QED system becomes 



H — J d 3 xip(x)(— «7 • V + mo)'ijj(x) 

+ J d 3 x^(x) 7 ^(x)(i M (x) + Af{x)) + H e 



(6) 
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where, A ext is a classical field, while others are quantized. H em is a Hamiltonian for electro- 
magnetic field. If we descibe the interaction between fermions by using a photon propagator, 
we get 



H' = J d 3 xii(x)(-tf • V + m + ^A ext (x))ip(x) 



2 

+y / d 3 x J dS^xh^D^x - y)kv)l v ^{y)- (?) 

Here, D^ v {x — y) is a full photon propagator. Quantized fermion- field operators ^p(x) and 
ip(x) can be expanded by one-particle state functions as 

Hx) = j2(4 + \x)d l + 4~\x)bt), (8) 

i 

kx) = J2(4 + \x)dl + 4- ) (x)h i ), (9) 

i 

where +(— ) means the electron ( positron ) state, % denotes quantum numbers for one- 
particle state, aj(dj) is electron creation ( annihilation ) operator, bf(Sj) is positron creation 
( annihilation ) operator, and tp^ is the four-component bispinor given like 

,(±) _ { A±) A±) ,(±) A±) )T , w) 
ri ~ h{2)i Ji(3)' Ji(4)^ ' ^ iu J 

where, T indicates the transposition of a matrix. 

After these preliminaries given above, we investigate the QED-DMFT. For our aim, first 
we introduce the one-particle local density matrix: 

Q(x) 4 x4 = -(^(x)^(x)) 4 x4- (11) 

Here if) and ijj are usual Dirac field, and they are 4-component bispinors. Thus our density 
matrix is 4x4 matrix, as denoted above. We can expand the 4x4 density matrix into the 
16-dimensional complete set of 7-matrices: 

16 

Q 4X 4 = £ V aQ A = Q S I + Qll" + Q^cT + C#7fcV + Q P ry 5 , (12) 

A=l 

Q A = - A tvT A Q = - A (i,T A iP). (13) 
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In this expansion, we take a convention as S denotes the scalar, V denotes the vector, T 
denotes the 2-rank antisymmetric tensor, A denotes the axial vector and P denotes the 
pseudoscalar. I is the 4x4 unit-matrix, f is usual Dirac gamma matrix, a^ v is defined as 
a^ u = \[f, r ) v ], and 75 = ?7°7 1 7 2 7 3 . Hence we obtain the Lorentz structure in our density- 
matrix Q. Each component of the expansion in (12) is given as follows: Q s = \trQ = \{tpip), 
Ql = \trfQ = \(4>f4>), Qlu = = i(V^V>>, Qi = \try 5 fQ = 1(^757^> and 

Q p = ^trij 5 Q = \{ipi^ip). Therefore, Q(x) includes a larger amount of informations than 
4-current j^x) = —e^^^ip^x)). It is clear from the discussion given above, if Q(x) is 
given, then and also other components can be calculated. 

The Rayleigh-Ritz variation principle of QED-DMFT should be given rather straightfor- 
ward extension of available theory, like the situation of the generalization of the Hohenberg- 
Kohn theorem to 4-current RDFT [7,14,18,25,26]. The exact energy functional of N- 
representable Q(x) for ground state of a QED system ( Strictly spearking, the particle 
number iV is not conserved in QED, but the total charge is conserved, as denoted in (15). 
) is given in the following form: 

E[Q(x)\ = (m Q \H'\m Q ) 

= F[Q(x)\ + J <Pxj»(x)Af(x). (14) 

Here, we use the Hamiltonian H' given in (7). Now A e * l (x.) is a static external field. Here 
we only consider a stationary problem. Variational search of the minimum point of the 
functional should be performed under constraint: 

dp f{x) = 0, j d 3 xj (x) = const. (15) 

First equation is arised from the gauge invariance of the theory derived by Nother theorem, 
while second one expresses the charge conservation. The fundamental theorems of DFT 
proved by Levy are constructed by two points [14] : (i) The variation principle can be applied 
to the ground state energy functional E[p(x)] or E[y(x,x')], where p(x) is the one-particle 
density while 7(2:, x') is the one-particle ( so called first-order ) density matrix, (ii) The 
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ground state of a system can be represented by p(x) or 7(2;, x') for ground state ( pgs( x ) or 
1gs{ x i x ') )• K is clear from his logic, the theorems do not distinguish between nonrelativistic 
case and relativistic case, and can be applied to relativistic case. Therefore, we argue that 
our functional E[Q(x)] will obey the variation principle and Q(x) for the ground state ( 
Qgs(x) ) represents the ground state of a QED system. Thus 

F[Q] = mm(^\f+V f ^}, (16) 
E[Q\ = min{F[Q] + f d 3 xf(x)A^x) J d 3 x Jo (x) = const.}, (17) 

where T is the kinetic energy while Vf-f is the interaction between fermions, presented in 
the Hamiltonian (7). 

Next step of our plan is the investigation into the Kohn-Sham scheme in our theory. The 
energy functional is devided into the following form: 

E[Q{x)\ = T s [Q{x)\ + J d 3 xf(x)Af(x) + \jd 3 xj rfV ^*^* + E xc [Q(x)}. (18) 

Here, T S [Q] is the one-particle kinetic energy functional 

T s [Q] = Jd 3 x Yl iHfV + mo)^, (19) 

-m<e n <e F 

and the 4-current j M (x) is given as 

j M (x) = -e Yl ^n(x)7^^„(x), (20) 

~m<e n <ep 

where, ip n is the Kohn-Sham one-body function in our theory. E X c[Q( x )] is the exchange- 
correlation energy functional in our theory. The definition of E X c[Q(x)} is given as 

E xc [Q{x)\ = F[Q{x)\ - T s [Q{x)\ d 3 x J d 3 x' f ^^\ (21) 

The third term of the right hand side of (18) is a Hartree-like term. This term arises if we 
replace D^ v (x — y) of (7) to the 0-th order photon propagator of the Feynman gauge: 

r dkn , s p*|fco||x-y| 

D '> ~ »> = / ^' tM ^k^7\ ' (22) 



Here g^ is usual metric tensor and defined as g^ = diag(l, —1,-1,-1). Therefore, the third 
term of (18) includes gauge parameter. E X c has to include contributions to compensate the 
gauge dependence of the Hartree-like term, for retaining the gauge invariance of the theory. 
From the variation with respect to the one-particle function, we obtain the Dirac-Kohn-Sham 
equation in our theory: 

16 

(-ry • V + m + ]T T A v A {x))^ n {x) = e nl °^ n (x), (23) 

<(x) ^ r(x) + /A ,^ + «l, (25) 

t ,s SE xc \Q(x)} 
V(X> = fp ' <26> 

<w ■ ^gp < 27 » 

Because our energy functional depends on the matrix Q(x), generally we obtain 16- 
component fictitious potentials. From the symmetry consideration, it is clear that, scalar 
potential v s gives the mass correction 5m, which will arises from radiative corrections, 
appears also in usual 4-current RDFT [7,26]. In fact, QED-DMFT is a non-perturbative 
method in QED. It is a famous fact that, QED becomes finite if its calculation fulfills the 
requirement of the gauge invariance [34] ( though we should employ gauge-invariant regu- 
larization and renormalization ). Therefore, if we can perform a calculation under retaining 
gauge invariance exactly, we will obtain a finite value for QED-DMFT. 



B. Relations between various DFT Schemes 

Now, we discuss the relation between various DFT, as depicted in Fig. 1. If we take 
into account only the components of 4-current in the expansion (12) of Q(x) for QED- 
DMFT energy functional, we obtain quantum electrodynamical current density functional 
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theory ( QED-CDFT ) [7,25,26]. Next, under performing the Gordon decomposition to 
spatial components of current j(x) = I + x S ( here I = ^{"^V^ — (V^)^} is the 
orbital current density while S = ^ - 0757°7'0 is the spin density ), and then completely omit 
I, we obtain relativistic spin density functional theory ( RSDFT ) [25,28]. Next, neglect 
the contribution of S for RSDFT, then we obtain relativistic density functional theory ( 
RDFT ). Performing some nonrelativistic reduction, especially only take into account the 
contribution of the large component to RDFT, and neglect all relativistic corrections ( spin- 
orbit, Darwin, mass- velocity etc. ), then we obtain nonrelativistic density functional theory ( 
NRDFT ) [5~7]. Generalizing the NRDFT to spin polarized system, we get nonrelativistic 
spin density functional theory ( NRSDFT ) [5~7]. Add the interaction between orbital 
current and external vector potential A ext to NRSDFT, we obtain nonrelativistic current 
density functional theory ( NRCDFT ) [35]. 



C. QED-TDDMFT 

Extension of QED-DMFT to the time-dependent region can also be considered. Starting 
point here is the action integral for density functional, like the QED-TDHF [31]: 

S[Q] = jf ' dt(^{t)\i^ - H'\**(t)). (29) 

Here matrix Q(x) = Q(x Q ,x.), depends on time xq. H' is given as (7). We take the 
Dirac-Frenkel variation principle to produce a stationary point at time- dependent density 
Q(x) [17,31]: 

SS[Q] = Sfj dt(*Q(t)\il - H'\* Q (t)) = 0. (30) 



Therefore we obtain the Euler equation: 

5S[Q] 



0, A=l~16. (31) 



5Qa(x ,x.) 

The Hohenberg-Kohn theorem in time-dependent potential was discussed by Runge- 
Gross [16,17] for nonrelativistic case, and it was extended by Rajagopal [30] for relativistic 
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case, in detail. In fact, we simply assume that our functional of Q can be used to discuss 
variation principle parallel with their results. 

The time- dependent Dirac-Kohn-Sham equation in our theory becomes 

16 



(i^dp -mo-Y, ^AV A {x))^ n {x) = 



|X — X 

T 



v p (x) = 



5QV(x) 
5E xc [Q(x)] 



SQT u (x) 
5E X ciQ(x)} 



SQ$(x) 
5E xc [Q(x)} 



(32) 

(33) 
(34) 
(35) 
(36) 
(37) 



SQ p {x) ' 

Here the effective single-particle Kohn-Sham potentials are also presented. Now, all of them 
are time-dependent potentials. 



III. GROUP THEORETICAL CLASSIFICATION OF QED-DMFT SOLUTIONS 



A. Group theoretical Classification 



In this section, we give a discussion of group theoretical classification of the QED-DMFT 
solutions. 

In the relativistic theory, we usually treat the Poincare group ( 4-translation and the 
Lorentz group ), charge conjugation, parity ( spatial inversion ) and time- reversal. Under 
the charge conjugation, ip and tfj are transformed as 



(38) 



where T denotes the transposition of a matrix, C = i^f 2 ^ is the charge conjugation matrix. 
Then Q(x) is transformed as 
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Q = -(W>) -> -C{^) T C~ l = CQ T C- 1 

= Q S I - ~ Ql^ v + Qtl^ + Q P H 5 - (39) 

Under the time reversal, 

^(t)-TV(-t), $(t) ^ $(-t)T. (40) 

Here T = 27 1 7 3 . Then, together with the rule of taking the complex conjugate about 
c-numbers, we obtain 

Q(t) = -(mm) - -T^{-t)mt))*T = TQ(-tyr 

= Q s *(-t)I + Q y *H) 7 ° - Q V *(~t) -l + Ql*(-t)a 0i 

-Qj*(-ty> + Q£*H)757 - Q A *(~t) ■ 757 - Q P *(-t)*75, (41) 

where i,j = 1,2,3. Qq is 0-th component of vector while Qq is 0-th component of axial 
vector. Q v and Q A are space components of vector and axial vector, respectively. Under 
the spatial inversion, 

V>(x) - 7 V(-x), ^(x) - i/S(-x)7°, (42) 
is transformed as 

Q(x) = -(^(x)V^(x)) - - 7 o^(-x)^(-x)) 7 ° = 7°Q(-x)7° 

= q 5 (-x)/ + g y (-x) 7 ° - q v (-x) ■ 7 - g^(-x)^ 

+Q5(-x)(7« - Q^(-x)757° + Q A (-x) • 757 - Q P (-x)*75- (43) 

The QED-DMFT solutions can be group theoretically classified into several types. To 
consider this problem, we determine the symmetry group of a system. In the atomic or 
molecular systems, the translation invariance is broken. In the case of an atom, only 0(3) 
rotation remains in the Lorentz group ( In the case of a molecule, 0(3) is replaced by 
the point group. ). Under the 0(3) rotational symmetry, we expand the Q s , Q p , Qq 
and Qq by the scalar spherical harmonics, while we expand the Q y , Q A and Q 1 - by the 
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vector spherical harmonics [36]. It is clear from (41) and (43), the behavior of each type of 
the density matrix given in (12) under the spatial inversion and time reversal depends not 
only on the structure of the 7-matrix, but also on the angular momentum of the spherical 
harmonics. Let us consider the case of an atom. We treat the group G = 0(3) x P xT. We 
introduce the subgroup of G as 

0(3) xPxT, 0(3) x P, 0(3) x T, 

PxT, 0(3), P, T, I. (44) 
For example, 0(3) x P x T-invariant solution is given as 

Q S * = Q S , Q%* = Q%, others = 0. (45) 
There is no case for 0(3) x P-invariant solution. 0(3) x T-invariant solution is given as 

Q S * = Q S , Qo* = Qo, Qo* = Qt others = 0. (46) 
0(3)-invariant solution is given as 

Q S * = Q S , Qo* = Qo, Qt = Qt Q P * = Q P , others = 0. (47) 

Due to the 0(3) rotational invariance, each density matrix can only have an s-wave compo- 
nent in all cases (45)~(47). All solutions given above are the cases for closed shell states. 
We may have magnetic QED-DMFT solutions for systems under discussions, if time reversal 
symmetry is broken ( for example, the case of (47) ). 

Under the presence of the vectorial density matrices like the Q v , Q A and QF V , the 0(3) 
rotational symmetry will be broken. For P x T-, P-, T- and 1- (no symmetry) invariant 
solutions, simple classification is impossible, because there are various possibilities of the 
angular momentum dependences of Q. This situation demands us futher investigation in 
more detail in the future. It must be noted that we have to solve the Dirac-Kohn-Sham 
equation to find which type of the density matrices ( and solutions ) will be realized. 



12 



B. Nonrelativistic and Ultrarelativistic Limits 



To consider the nonrelativistic and ultrarelativistic limits, we take the standard repre- 
sentation. The four component bispinor if) is partitioned as [36] 

(A 

1> = , (48) 

w 

where is the large component and x is the small component. Then we obtain 

Q s = ^trQ = i(^> = " X f X>, (49) 

Qo^ = \trfQ = I (^7 V> = + X f X> , (50) 

Q y = ±tr?Q = \{^) = \( X ^4> + ^x), (51) 

Qj v) = hra^Q = ^7°70> = i(-X^<P + f <?X>, (52) 



Qf A) = -tra^'Q = i (^757°7V') = - X^x), (53) 

1 n„ 1,t n 1 



= i tr 75 7°Q = i (^757» = - i (xV + ^ f X>, (54) 
Q A = ^tr 75 7Q = 1 (^75#> = - ^ (</> f ^> + X^X) , (55) 
Q p = ^tri 75 Q = ^(^75^) = \{~iX ] <P + itfx)- (56) 

Here, denotes vector-like components of Q? v , while denotes axial- vector-like com- 
ponents of Q T . Therefore, at the nonrelativistic limit ( x wm vanish ), the Q s and Qq 
coincide with each other and they give while the QJa) gives and the Q A gives 

— \$o<$>. The Q y , Q(y), Qq and Q p vanish at the nonrelativistic limit. (ft<f) corresponds 
to the usual density of nonrelativistic theory, while <$(j§ corresponds to the spin density of 
nonrelativistic theory. Futhermore, at the ultrarelativistic limit ( <fi and x coincide with each 
other ), Q%, Q y , Qq and Q A remain while others vanish. Therefore, for a system under 
strong relativistic effects, Q y and Q A may give important contributions to QED-DMFT 
solutions. 
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IV. DISCUSSION 



In this paper, we have investigated the QED-DMFT and the group theoretical classifi- 
cation of its solutions. We have introduced the one-particle local density matrix Q(x), and 
have decomposed by using 7-matrices. We have discussed the transformation properties of 
the matrix under group operations, and its nonrelativistic limit. 

Now, we discuss the application and extension of our theory. 

In the context of the electronic structure of atoms, there are three effects: The electron 
correlation, the relativistic effect and the QED effect. The electronic structure of atoms is 
determined by the relation of these three factors. The electron correlation depends on the 
electron numbers and can be treated by several many-body techniques. The relativistic effect 
( kinematic effect, spin-orbit, Darwin, mass- velocity, etc. ) becomes large with increasing 
the atomic number. On the other hand, the case to describe the inner core electrons of 
heavy elements or, electrons of highly ionized heavy elements such as lithium-like uranium, 
the QED effect ( mainly originate in negative energy state ) can not be neglected and we must 
take the Dirac sea into account. In the case of heavy atoms, as the ionicity becomes high 
and the electron number decrease, the many-body effect becomes small and the QED effect 
becomes large. In principle, QED-DMFT can exactly estimate the three effects discussed 
above. Therefore we propose that QED-DMFT should be applied to the cases where both 
the many-body effects and the QED effects can not be neglected. Practically, we need a 
well-behaved exchange-correlation functional which contains both the correlation energy and 
QED effect satisfactorily. Heavy elements in middle level of ionicity should be one of the 
subject for our theory. The calculation of electromagnetic properties, for example, g-factors, 
hyperfine interactions, nuclear magnetic resonance shielding constants, are interesting. The 
collision of two uranium atoms is also one of the interesting subject of QED-TDDMFT. QED- 
TDDMFT should be useful to describe excitations, susceptibilities and collective modes of 
systems. Recently, experiments of x-ray irradiation to cluster plasma are performed, and 
verious new phenomena were studied [37]. The cluster of heavy ions under middle revel of 
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ionicity is now an interesting subject in this area. The importance of the relativistic and 
QED effects is discussed in such objects. 

It is a famous fact that, the available exchange-correlation functionals only describe 
the dynamic correlations reliably, and they fails for systems with strong static correlations. 
Thus, usual Kohn-Sham scheme often breaks down in the case of near degenerate ground 
states. In the usual method of quantum chemistry, the static correlation should be described 
by using multi-determinantal wavef unctions. The CASSCF is one of the methods to take 
into account the static correlation effect correctly. Recently, the CAS-DFT method was 
proposed [38]. In this method, the dynamic correlation is considered by usual DFT scheme, 
while static correlation is handled by CAS-type wavefunction. The energy functional is given 
as 

E[p] = min{F[p] + J d 3 xp(x)v(x)} , (57) 

F\p] = F2f\A + E% AS \p], (58) 
F%f[p]= min (* CAS \f + V ee \V CAS ), (59) 

where T and V ee are the kinetic energy and electron-electron interaction, respectively. ty CAS 
is a CAS trial wavefunction. In the scheme given above, exchange energy and static cor- 
relation energy are covered by F£ AS , while E^, AS includes the dynamic correlation energy 
only. CAS-DFT in nonrelativistic case gives accurate results [38]. Thus we propose that 
this method can be extended to the relativistic theory, and should become usuful tool to 
consider near degeneracy effect in relativistic DFT framework. For this aim, we introduce 
the density matrix which is defined as 



^ ++ H — \ 
Pij Pij 



^ ($os|a]ai| $ Gs) {$Gs\bjai\&Gs) ^ 



\ Pij Pij J 



(60) 



y ($Gs\a]b\\$Gs) ($Gs\bjH\$Gs) ) 

Then we diagonalize the density matrix, we obtain the natural orbital given by four- 
component bispinor rff 1 ^ and 7ij as its occupation number. By using the occupation number 
rii, we can select the active space for treatments of the near degeneracy effects [39]. The 
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natural orbital analysis of a resulted solution is useful to elucidate the type and magnitude 
of broken symmetry. 

It is also interesting for us to use the information entropy to measure the electron corre- 
lation in various systems [40]. For example, Shannon entropy 

S Shannon (p) = - J d 3 xp(x) In p(x) , (61) 

or Jaynes entropy 

S'*"™ 8 = -'EmhMi, (62) 

i 

might be useful tools. The electron correlation energy can be estimated based on the Collins's 
conjecture: 

E corr = -kJ2ni In m oc S Jaynes . (63) 

i 

Therefore, by using these entropies, we propose to estimate the strength of the electron cor- 
relation in various atoms under various ionicities with QED-DMFT. The relations between 
various theories proposed in this paper are depicted in Fig. 2. 
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FIGURES 

FIG. 1. The relations between various DFT. Meanings of abbreviations are given in text. 
FIG. 2. The relations between various theories proposed in text. 
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